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I . Half a century ago, Casimir— predicted that the quantum fluctuations of the electromagnetic field within a planar 
] cavity would produce an attractive macroscopic force on its boundaries. His prediction was based on the properties of 
the field when confined by perfectly reflecting mirrors. It is only recently that experimental studies have attained the 
necessary accuracy to test in detail the theoretical predictions;^ the Casimir effect has now been measured'^'^"^'^ with 
uncertainties as small as 1% and at distances^ down to « 60nm. Therefore, theories of the Casimir force that account 
for the properties of realistic cavities have become indispensable. In 1956 Lifshitz proposed a macroscopic theory for 
two semi-infinite homogeneous dielectric slabs^ characterized by complex frequency dependent dielectric functions. 
. The stress tensor was obtained from the self-correlation of the fluctuating electromagnetic field whose source consists 
^ ■ of fluctuating charge and current densities within each slab. Their autocorrelations are related through Kubo's 
O^i formalism,^ causality, and the fluctuation-dissipation theorem, to the complex dielectric response. In his calculation, 
Lifshitz considered only homogeneous and isotropic media and it was assumed that fluctuating sources at a given 
, position were completely uncorrelated with sources at nearby positions so that the results were directly applicable to 
J> ' semiinfinite homogeneous flat local media, and did not cope with more complex systems, such as thin films, layered 
OO ' systems, superlattices, photonic crystals and metamaterialsJ^ For the same reason, it seemed incapable of dealing 
■ with spatial dispersion and screening at realistic surfaces ^iiii^ii^ 

' Several alternatives to the derivation of Lifshitz have been proposed. Barash and Ginzburg^^ determined the 
, allowed frequencies of the cavity modes by solving Maxwell's homogeneous equations and imposing planar boundary 
0^ conditions. The energy could be obtained from the resulting density of states if there were no dissipation. Nevertheless, 
, dissipation yields complex frequencies and the interpretation of the solutions as normal modes looses meaning. Barash 
and Ginzburg overcame this problem by introducing an auxiliary non-dissipative system. The use of an auxiliary 
system was further developed by KupiszewskaJ^ in a ID calculation in which the problem of quantizing a dissipative 
^ ' system is attacked by accounting both for the dynamics of the vacuum modes and of the atomic dipoles to which they 
k>( \ couple and which make up the material, together with a thermal reservoir in which the atomic radiators dissipate 
. the absorbed energy. A disadvantage of this approach is that it requires an explicit microscopic model for the walls 
d ' of the cavity and for the thermal bath, thus appearing to restrict its generality. A similar approach™ was based 
on a Green's function method and Kubo's theorem. In both cases the stress tensor is obtained from the vacuum 
modes with an explicit dependence on the dielectric response e(a^). An alternative treatment of the Casimir force 
was introduced by Jaekel and Reynaud^^, who calculated the radiation pressure within a cavity bordered by partially 
transmitting but lossless mirrors. Each mirror was replaced by an infinitesimally thin scatterer characterized by a 
unitary, energy conserving scattering matrix. Their calculation was later generalized to the case of lossy optical 
cavitieai^ by complementing the cavity modes with noise modes in such a way that the total scattering matrix was 
unitary. The scattering matrix corresponding to the cavity modes was then obtained through the optical theorem. 

Mochan et alr^t^ have obtained an expression for the Casimir force using both the scattering approach and a 
dissipationless ancillary system. They have argued that in thermal equilibrium, all of the properties of the radiation 
field within a cavity are completely determined by the optical reflection amplitudes of the walls. Thus, the Casimir 
force may be obtained from the stress tensor of any system whose reflection amplitudes are identical to those of 
the real system. A dissipationless fictitious system with those properties was conceived: It had infinitely thin walls 
characterized by a unitary scattering matrix whose elements corresponding to the optical reflection amplitudes from 
within the cavity were chosen to be identical to those of the real system. The transmission amplitudes were chosen in 
such a way that the energy that was not reflected was transmitted without loss to the vacuum outside of the cavity. 
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This permitted a full quantum mechanical calculation of the fields, even when the real system is dissipative. The 
field modes were quantized and counted by adding perfect mirrors far away from the walls of the real cavity. These 
quantizing mirrors produce a field that mimics the incoherent radiation back into the cavity that is responsible for 
maintaining a detailed balance and thus the thermodynamic equilibrium in the case of lossy or dissipative real mirrors. 
The field that enters the cavity after being reflected by the quantizing mirrors has a very large, frequency dependent 
phase that becomes infinitely large as the mirrors are moved infinitely far away. 

The main result from the work mentioned above is that if Lifshitz formula is written in terms of the reflection 
coefficients of the walls of the cavity, or equivalently, in terms of their exact surface impedance2ii2^, it becomes 
applicable to any system with translational invariance along the surfaces and isotropy around their normal and not 
only to semiinfinite, homogeneous, local mirrors. Thus, it may be employed to calculate the Casimir force between 
semiinfinite or finite, homogeneous or layered, local or spatially dispersive, transparent or opaque, finite or semi- 
infinite systems. Through a simple substitution of the appropriate optical coefhcients, the formalism has allowed 
the calculation of the Casimir force between photonic structures ^^24 non-local excitonic semiconductors non-local 
plasmon-supporting metals with sharp boundaries,— and between realistic spatially dispersive metals with a smooth 
self-consistent electronic density profile. The relative simplicity of the formalism has allowed its generalization to 
non-isotropic systems and the calculation of Casimir torques. With a few modifications, it has also been employed 
for the calculation of other macroscopic forces, such as those due to electronic tunneling across an insulating gap 
separating two conductors. 

Nevertheless, there was a shortcoming in the derivation of the Casimir force presented in Refs. T^l20j as it was 



uncritically assumed that a unitary, energy-conserving scattering matrix could be built through a proper choice of 
transmission coefficients. Somewhat surprisingly, it turns out to be impossible to find a unitary scattering matrix 
for evanescent waves, i.e., for Q > w/c, where Q is projection of the wavevector parallel to the cavity walls, oj is 
the frequency and c the speed of light, as a single transmitted evanescent wave is unable to transport energy away 
from the surface while an incident and a reflected evanescent wave do transport energy from the cavity towards the 
surface of lossy and dissipative systems.— Although the contribution of evanescent waves to the Casimir force could 
be obtained as an analytic continuation from the region of propagating waves, it is not obvious a priori that this 
extrapolation would yield the correct result. 

Recently?^ the problem of energy transport in the evanescent region was dealt with by modifying the fictitious 
system introduced in Refs. [T9ll23| in such a way that evanescent waves within the cavity couple to propagating waves 
outside the cavity. This was accomplished by filling completely the region beyond the infinitcsimally thin mirrors 
with a dispersionless and dissipationless fictitious dielectric with a large permittivity ey. In this way, the fictitious 
light cone Q < ^/ej:Lu/c extends beyond the light cone Q < w/c of the vacuum cavity and in the limit e/ ^ cx) all 
of the evanescent waves within the cavity would be able to couple to propagating waves within the fictitious region. 
Lifshitz formula was thus proven to be valid both within and without the light conci^ 

One drawback of the calculation presented in[3ll is its use of an extremely unrealistic dielectric, with a suspiciously 
large, real, frequency independent dielectric constant e — > oo. As the Casimir force ought to be determined by the 
refiection coefficients of the real mirrorS )^^'^° all of the details of the fictitious system beyond the walls of the cavity 
ought to be superfluous after their contribution to the reflection amplitude has been accounted for, and it should be 
possible to setup the calculation without the need of specifying them. 

The purpose of the present paper is to develop yet another derivation of the Casimir force within cavities with 
walls made up of arbitrary materials characterized only by their optical reflection amplitudes. As in Refs. ^19,^20 , jilj] , 
we introduce a dissipationless fictitious system with no degrees of freedom beyond those of the electromagnetic field 
and with a cavity whose walls have the same optical coefficients as the real system. However, unlike the calculations 
above, we avoid giving any detail of the fictitious system beyond the reasonable fact that it should be consistent 
with detailed balance so that thermodynamic equilibrium is satisfied, i.e., on the average, for each photon that is not 
coherently reflected at a cavity wall and is therefore either absorbed or transmitted beyond the system, an identical 
photon has to be incoherently injected back into the cavity with no phase relation to the lost photon. We believe that 
this derivation of the Casimir force is quite simple and that it can be readily generalized to other geometries, allowing 
the calculation of the dispersion forces in cavities of varied shapes whose walls are made up of realistic materials. 

The structure of the paper is the following: First we review briefly the model employed in Ref. [Slj . In Sec. |TI] we 
study the scattering matrix of a fictitious interface between vacuum and a dispersionless dielectric with an additional 
infinitcsimally thin scatterer that forces the optical coefficients to be the same as those of the real system and in Sec. 
mil we study the reflection amplitudes for both propagating and evanescent waves after adding quantizing mirrors. 
Then, in Sec. IIVI we eliminate the superfluous details from the calculation, keeping only the delay T oo before 
injecting back into the cavity the photons it looses through absorption or transmission in order to restore equilibrium. 
In Sec. |V]we obtain the electromagnetic normal modes of the cavity and the contributions of the cavity walls to the 
density of states which we employ in Sec. IVII to obtain their contribution to the thermodynamic properties. Finally, 
section Ivn] is devoted to conclusions. 



3 



Vacuum: e = 1 / Dielectric: £ = £/>! 




FIG. 1: Infinitely thin scattering interface (/) between vacuum and a dispersionless dielectric with permittivity e/. Incoming 
(j) and outgoing (o) waves within vacuum (v) and the dielectric (d) are illustrated, as well as the projection of their wavevectors 
parallel (Q) and normal (A:„ and kd) to the surface. The choice of coordinate axes is also shown. 



II. ENERGY FLOW 



In Ref. ^31, a fictitious system was introduced, consisting of a vacuum cavity bordered by two infinitesimally thin 
sheets followed by dispersionless dielectric slabs which are terminated by perfectly reflecting mirrors. It was argued 
that in equilibrium the electromagnetic field within the real cavity ought to coincide with the electromagnetic field 
within the fictitious cavity. The reflection amplitude of the infinitesimal sheet together with the dielectric was chosen 
to coincide with the reflection amplitude of the mirrors that make up the real cavity. The perfect mirrors were 
incorporated in order to quantize the normal modes and in order to inject back into the cavity any radiation that is 
not coherently reflected at the surface, guaranteeing thermodynamic equilibrium. The re-injected radiation acquires 
a large phase as it travels across the wide dielectric slab, thus mimicking the incoherent re-radiation of photons that 
are lost through transmission or through absorption at the walls of the real cavity. 

Consider one of the mirrors of that flctitious cavity, as illustrated in Fig. [2 consisting of an infinitely thin reflector 
between vacuum and a dispersionless dielectric. We first assume that all the waves in the figure are propagating, i.e., 
ky and are real, then the energy flux along the normal, taken to be the z direction, is 

2 2 

Sy, = ^fc.(N.|' - \oX), Sdz = T^fcd(|od|2 - \idf) (1) 

OTTCJ OTTCJ 

in the vacuum and the dielectric side respectively, where ia and Oq are the amplitudes of the incoming and outgoing 
electric fields within the a — v,d region. Energy conservation requires that 

kv\ov\^ + kd\od\^ = ky\iy\^ + kd\id\^ (2) 
for arbitrary values of the incoming amplitudes iy and id- If we define a scattering matrix S through 

/ VkyOy \ ^ a f ^ 



then S should be unitary in the usual sense, i.e., S^S — SS^ — 1, where 1 is the unit matrix. We identify the 
components of S as 



rv td.y'ky/kd 

tv^/kd/ky rd 



(4) 



where and ta are the reflection and transmission amplitudes corresponding to incidence on the interface from 
medium a. Unitarity then yields the relations 



|2 I U |2 _ 1 I |2 I I . |2 _ 1 * . / ^t' I .* r^d 



-^\ty\'^l, \rd\' + ^\td\' = l^ r:tdJ^+rdt:J^=0, (5) 

Ky Kd V V >^v 



which imply 



Rv = Rd, Ty = Td, Ry + Ty — Rd + Td — 1, (6) 



e = l / £ = £/ M 
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FIG. 2: Multiple reflections at a fictitious mirror made up of the interface shown in Fig. [T]with the dielectric terminated by a 
perfect mirror (M) after a distance Ld- 



where we identify as usual the reflectance — Rd = kdP) and transmittance [Ty — (fcc;/fc«)|i«P, Td — 

{ky/kd)\td\^) and we denote by (. . .)* the complex conjugate of any quantity (...) . Of course, Eqs. ([6]) are consistent 
with Fresnel relations^. However, we must remark that and ta above are not given by the Fresnel relations due 
to the presence of an additional scatterer at the interface which forces the reflection amplitude from the vacuum side 
to agree with the reflection amplitude of a boundary of the real cavityii2i2£ii^ 

In the case of evanescent waves, i.e., when Q > uj/c the normal component of the wavevector fc„ — iKy is imaginary 
and the energy flux is not given by Eq. ([1]) but by 

^2 

KyiilOy)". (7) 



Sttui 



instead of Eq. ([T]), where we denote by (. . .)' and (. . .)", or equivalently, by Re(. . .) and Im(. . .) the real and imaginary 
parts of any quantity (. . .). If Q were so large Q > ^JtjLojc so that waves were also evanescent in the dielectric, we 
would also have 

Sd. = 2—Kd{idO*d)" (8) 

and it would be impossible to satisfy energy conservation S^z = Sdz for arbitrary incoming ainplitudes. This was 
the reason for introducing the fictitious dielectric e/ and for taking the limit — > oo in Ref. |3l| |. For waves that 
are evanescent within vacuum but are propagating within the dielectric, Sdz would be given by Eq. ([1]) and energy 
conservation could be satisfied provided 

|r,|2 = l, |i„|2^2^r;', U = ^^tyd. (9) 



III. TOTAL REFLECTION AMPLITUDES 

If we terminate the dielectric in Fig. [T] a large distance Ld from the interface and situate a perfect mirror there, 
then the total reflection amplitude of the system may be obtained by adding multiple reflections (see Fig. [2]), 

n^ry- Utd^^^-^-{l - rde^^^-^- + r^e^^^^^- -...) = + ~ ^iM"""''^' ■ 

1 + TdC'^^'^'i^'i 

If the waves in vacuum are propagating, we may use Eqs. © to obtain 

n = , ■ (11) 

1 + Tde 

Although the reflection amplitude of the real system |r^| is typically smaller than 1, it is easy to verify employing 
Eqs. (O that |rtp = 1, so that all of the energy that crosses the interface / eventually crosses back after some delay. 
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On the other hand, for evanescent waves we can use Eqs. ([9]) to write 

so that although r„ is complex, rj is a real quantity. According to Eq. ([H]), this means there is no net energy flux 
towards the interface so that energy balance would be achieved within the vacuum cavity also for evanescent waves. 

As argued in Refs. [19,20,31], the properties of the electromagnetic field within a cavity made up of lossless mirrors 
such as those described in this section ought to agree with those within the real cavity. Thus, equations equivalent to 
pT|) and (fT2)) were employed in Ref. [3l[ to calculate the normal modes of the cavity and from them its thermodynamic 
properties, including the Casimir force. The purpose of this paper is to discard the unnecessary and dubious details 
of the fictitious system, such as the dispersionless dielectric with a large permittivity. 

IV. DELAY 

Eqs. pT|) and (fT^ contain elements from the real system, such as the reflection amplitude r„ of the cavity walls. 
They also contain quantities that relate to fictitious quantities, such as the dielectric constant £/ of the dielectric, 
its width Ld and the reflection amplitude for light impinging on the interface from the dielectric. To eliminate 
those quantities from the model we first notice that Eqs. pT|) and contain terms which are relatively slowly 
varying functions of the frequency, such as the refiection amplitudes of the real cavity walls. On the other hand, 
they contain extremely fast varying functions of the frequency such as e^*'^'*^'', which are due to the long time taken 
by the field to transverse twice the width of the L4 dielectric after entering it across the interface in order to return 
after reflection in the perfect mirror. This long delay is precisely what allows the fictitious system to replenish those 
photons that are lost from the cavity with such a large phase that they mimic the incoherent thermal photons that 
would be radiated back into the cavity in the real lossy system. It is actually the essence of the fictitious cavity 
model: In the real cavity photons that are not coherently reflected are lost through transmission or absorption at the 
walls; the cavity walls would radiate photons incoherently and inject them into the cavity to sustain thermodynamic 
equilibrium. In the fictitious system photons are not lost, but they are delayed a time T — 2kdLd/i-u — > 00 before they 
reach the cavity again, so they are essentially indistinguishable from incoherently radiated thermal photons. We can 
keep the delay T in the model eliminating all other details by postulating that the reflection amplitude takes the form 

where a and b are slowly varying functions of the frequency to be determined. The term r„ in the numerator accounts 
for the coherent refiection of the real system. The term ae*""^ corresponds to re-radiation by the cavity walls or to 
thermal photons entering the system from outside to replenish the cavity losses, and it includes the long delay T 
which mimics incoherence. Finally, it may happen that such a photon is reflected back from the interface into the 
walls and doesn't enter the cavity on its first attempt. Thus, it could be absorbed and re-emitted some time later. 
The possibility of multiple re-injection attempts is accounted for by the term 6e"^"^ in the denominator, analogous to 
the denominator in Eq. (jlOp . 

In equilibrium all the energy that leaves the cavity has to enter it again. Thus, for propagating waves the total 
reflection amplitude must obey 

\rt\' = l, (14) 

which yields 

\r, f + \af + 2Re(r:ae'"^) = 1 + |6|2 + 2Re(6e'"^) (15) 



after substitution of Eq. (JTSj) . Separating the slowly and rapidly varying functions of lu in the previous equation we 
obtain the two equations 

6 = r>, |r,|2 + |a|2^1 + |6|2 (16) 

from which we obtain 



(17) 



z\ L zi 



FIG. 3: Cavity made up of two plane, parallel, lossy mirrors (1 and 2) with reflection amplitudes r\ and r-i situated a distance 
L apart at z\ and Z2. As the field is reflected back and forth (arrows) its amplitude diminishes. 



where the real phase b = 5{u}) is some slowly varying function of uo. Substituting Eqs. (fT7|) in (fTS]) we finally obtain 

^* = l + ^*e'(*+-^)' ^^^^ 
For evanescent waves, energy conservation requires 

r'l = (19) 

instead of Eq. (|14p . which implies 

Im[r,„ + ab* + (a - rlb)e''^'^] = 0. (20) 
As done above for Eq. (jlSp we separate fast and slow varying functions of u to obtain 

a = r*b, lm{ry + ab*) = 0, (21) 

from which wc obtain 

b = e'^, a = r*y\ (22) 
where b is again a slowly varying function of the frequency. Substituting Eqs. (|22p in (|13p we finally obtain 



_ + r:e'(^+"^) ^ ^ Re(r. + r;e'(^+"^)) 
^* 1 + e*('5+'^'^) |i + e*('5+'^'^)p ■ 

Notice that Eqs. (fTT|) and p2|) have the form of Eqs. (fT8|) and ([23|) respectively, where we identify e** = ^u/'^d for 
propagating waves and e"^ = for evanescent waves, as = 1 in the former case according to ([S]) and jr^l = 1 

in the latter case according to ([9]). However, we didn't have to postulate any extraneous fictitious dielectric slab nor 
the reflection amplitude rd from within it. We will take the limit T ^ oo and we will show that the unknown phase 
S is irrelevant. 



V. NORMAL MODES 



Given a real dispersive and dissipative system with a coherent reflection amplitude r„, in the previous section we 
obtained an expression for the total reflection amplitude rt that accounts for the coherent reflection and mimics the 
incoherent reflection in thermodynamic equilibrium through terms that oscillate rapidly as a function of frequency. 
Consider now a planar cavity bordered by two arbitrary material slabs with reflection amplitudes ri and r2 (Fig. 
[31). For a given Q and oj the field within the cavity would be proportional to ArC**^"^ + A;e~''^"^, with and 
Ai the amplitudes of the right and left moving contributions to the field. Applying boundary conditions at zi 
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(Are*'^"^i = riAie and at Z2 {Aie ^^^^'^ = r2^,e ^'^f^s) would yield the usual condition for the normal modes 

of the cavity, 



D==l~ rirae^**^"^ = 0. (24) 



Due to the losses at the mirrors, we expect that after externally exciting any of the modes obtained from Eq. 
it would decay in time (Fig. [3]). Thus, for a given real Q, Eq. p4)) would yield a set of complex frequencies whose 
imaginary part is related to the finite lifetime of each mode. Nevertheless, at thermodynamic equilibrium, the energy 
of each mode would be replenished through re-radiation, so that the equilibrium modes would be stationary instead of 
decaying. We may obtain the stationary modes by replacing the reflection amplitudes ri and r2 by the corresponding 
total reflection amplitudes ru and r^t obtained through equations such as p8)) and (p3|) . 

A = 1 - r-itrate^^'^-^ = 0. (25) 

As Tit and r2t are lossless, the frequencies lo(^ that solve Eq. ([^5]) for any real value of Q are necessarily real. Thus, 
we can straightforwardly apply quantum mechanical methods to the modes derived from (|25p . In particular, the 
contribution of the ^-th mode to the energy of the system is simply — {{ni) + l/2)huji, where (ne) is the average 
occupation number of the mode. 

For propagating waves |rif | — \r2t\ — 1 (Eq. HH)), so that Eq. (^5)1 may be recast as 

arg(ritr2te2*'="^) = 2tt£ (26) 

with £ an integer. Rewriting (fT8|) as 

where we introduce the subindex z = 1, 2 to denote each of the two mirrors, Eq. (I26p becomes 

2 arg(l + rie"*(*i+"^i)) + 2 arg(l + r2e"*('^^+"^^)) + luT + 5 + 2fc,,L = 2n£, (28) 

where T = Ti + T2 and S = Si + 62- Notice that as lo varies, l + rie^'(''i+"^i) moves counterclockwise in the complex 
plane around a circle centered at 1 whose radius |ri| < 1. Thus, it does not encircle the origin and its contribution to 
the phase in (j28p is bounded. The same happens with the second term in (|28p . Thus, the eigenfrequencies uji never 
get far away from the eigenfrequencies 

com = (27r^ - 2fc„L - S)/T (29) 

corresponding to vacuum, for which ri = r2 — 0. Notice that as T — > 00 successive frequencies approach each other 
and the density of states diverges. 

Consider now a small frequency range fl of size Auj centered at a given frequency ZJ, where Aw is much smaller 
than any characteristic frequency of the system. Nevertheless, as we will take the limit Ti,T2 00, we may assume 
that AujT ^ 1 so that the number N{il) of normal modes uJi within fl is large N{il.) ^ 1. Using Cauchy's argument 
principle, we obtain 

where 7 is a clockwise closed path that encircles and 

f{Lu) = (1 + fje'(''i+'^^i))(l + f*e'('52 _ _^ e*('5i+-'^i))(r2 + e»('52 W2))g2^fc„L ^3^^ 

has within 7 the same zeroes as Dt (Eq. (|25p ) and no poles. Here we introduced linearized and therefore analytical 
approximations f^, r2 to rl, r2, taking advantage of the smallness of Aw, so that / is an analytical function even if 
Dt is not. Choosing a as a path that goes from ZJ — Allj/2 to ZU + Aw/2 a small distance rj ^ below the real axis 
and returns a distance 77 above the real axis, we can rewrite Eq. (j30p as 

N{n) = -i^ r^^"'^ dLo log[/(w - zr;)//(w + Z77)]. (32) 
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Approaching the hmit Ti, T2 ^ oo before taking the hmit 77 — > 0, we may assume that jyTi, 77T2 ^ 1 even if 77 <C Aw, 
so that /(w — irj) ^^{s+uJT)+r|T ^~*~* _ g2ifc„L-j _|_ jjy) ^ 1 _ 7-ir2e^*'^"^. Subtracting the number of modes 

-/Vo(n) corresponding to vacuum, which may be obtained from Eq. (|32p by replacing ri and r2 by zero, we obtain the 
contribution /^N{^) of the cavity walls to the number of modes 

AN{n)^--f^ ^ dujlm^\og{l~rir2e^'''-^). (33) 

I" Ju~Au/2 duj 

Notice that all the terms in Eq. ([55]) arc slowly varying and, as we chose a very small Aw, the integral becomes trivial. 
Dividing the result by Aw we obtain the contribution of the walls to the density of states 

p= --Im^log(l -nrae^^'^-^). (34) 

TT duj 



A similar procedure may be employed for evanescent waves by substituting Eq. ((23)l instead of (|18p in ((25l) . The 
number of modes within the frequency range Q is again given by Eq. ()30|) . but choosing 



/(w) = (1 + e'('5i+-'^i))(i + e«(<52+c^T.)^ _ ^ f*e'(''i+"^i))(r2 + f;e'('52+-'T2)^g-2«„L^ ^35^ 

where we took = in^. Now we substitute /{cu-iri) gi(s+uiT)+,iT ^^-^^ _ ~*~*^-2k^l^ f{u! + ir]) 1 - rir2e'^''"-^ 
in Eq. pop and subtract from the resulting expression the number of modes corresponding to vacuum. The result is 
again given by Eq. ([55)1 , so that the contribution of the walls to the density of states for evanescent waves is also given 
by Eq. (p4)) . which is therefore valid both in the propagating (w > Q/c) and in the evanescent (w < Q/c) sectors. 

Finally, from Eq. (|34p we can obtain an expression for the contribution of the walls to the average of any real 
quantity 14^^(Q,w), namely 



(36) 



where 

Cm = (^^ar^2e^^'=-^)-\ (37) 

we incorporated the fact that light has two independent polarizations fi = s,p over which we summed, and we 
performed the usual sum over parallel wavevectors J2q • ■ • — * -4/(47r^) / cPQ . . . where A is the area of the mirrors. 

VI. THERMODYNAMIC QUANTITIES 

As simple applications of Eq. p6|) we calculate the ground state energy 

Uo = (?i^/2) = -^Im^ — Jd^Ql dww— log(l - C^'). (38) 

Performing the angular part of the wavevector integration and integrating by parts over frequency we obtain 

Uo = ^Ini^ dQ Q duj log[(l - C')(l - C^')]- (39) 



/o ^0 

Similarly, the internal energy at finite temperature 



U = (hiog) = Alm^ dQ Q dco —{cog) log[(l - C')(l " Cp') 



(40) 



where 



1 1 , f(3huj\ 1 
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is the occupation number of a state with frequency lo at temperature l/ksP, with fc^ Boltzmann's constant. The free 
energy 



(42) 



-^eschf^l (43) 



(44) 



where 

2 V 2 

is the partition function of a single mode of frequency lo. From (j42p and (|40p we may obtain the entropy 

^ = A^lmJ^^ dQ Q log[(l - C')(l - 

Finally, deriving the free energy with respect to L we may obtain the Casimir force 

which agrees with Lifshitz formula when written in terms of the reflection amplitudes. For the actual evaluation of 
the frequency integrals in Eqs. ([55|) - (|^5)) the integration path can be conveniently deformed from the real into the 
imaginary axis yielding the usual Matsubara summations. 



VII. CONCLUSIONS 



We have obtained expressions for the thermodynamic properties of a cavity formed by two flat mirrors. We 
derived our expressions without making any assumption about the mirrors, which were completely characterized by 
their reflection amplitudes ri and The mirrors could have been conducting or dielectric, opaque or transparent, 
dispersionless or dispersive, lossless or dissipative, semiinfinite or finite, homogeneous or layered, with abrupt or 
smooth boundaries, local or spatially dispersive, etc. To derive our results we postulated a fictitious system which 
has the same reflection amplitudes as the real system, but such that any energy that leaves the cavity eventually 
comes back but after a very long delay which is taken to infinity. Thus, we mimic in a closed system the incoherent 
field radiated by the walls of the cavity and the thermal radiation of the environment which replenishes the energy 
lost by the cavity in thermal equilibrium. The requirement of equilibrium allowed us to find total, lossless reflection 
amplitudes for the fictitious system which allowed a full quantum mechanical description of the system, from which 
we identified and counted the normal modes, obtaining an expression for the contribution of the walls of the cavity 
to the density of states of the system. With the density of states, obtaining expressions for all of the thermodynamic 
quantities becomes a straightforward task. The expressions we obtained agree with those found in the literature, 
although our derivation shows they are more general than implied by most other derivations. We believe that the 
procedure is simple enough to be easily extended to non-planar cavities made up of arbitrary materials. 
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